NX=0l & NY=0l

OPENR,LUN,/GET_LUN,'.\COEFF_DAY.DAT',/F77_UNFORMATTED
READU,LUN,NX
READU,LUN,NY
X=fltARR(NX) & Y=fltARR(NY)
READU,LUN,X,Y
BDXX=fltARR(NX,NY) & BDYY=BDXX & BDXY=BDXX
READU,LUN,BDXX
READU,LUN,BDYY
READU,LUN,BDXY
FREE_LUN,LUN

BDYY=BDYY/REBIN(REFORM(Y^2,1,NY), NX, NY)
BDXY=BDXY/REBIN(REFORM(Y,1,NY), NX, NY)

ALPHA=X*!radeg
EK=(SQRT(Y*Y+1D0)-1D0)*0.511
end